Phase-fluctuating 3D condensates in elongated traps 
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We find that in very elongated 3D trapped Bose gases, even at temperatures far below the BEC 
transition temperature T c , the equilibrium state will be a 3D condensate with fluctuating phase 
(quasicondensate) . At sufficiently low temperatures the phase fluctuations are suppressed and the 
quasicondensate turns into a true condensate. The presence of the phase fluctuations allows for 
extending thermometry of Bose-condensed gases well below those established in current experiments. 
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Phase coherence properties are among the most inter- 
esting aspects of Bose-condensed gases. Since the dis- 
covery of Bose-Einstein condensation (BEC) in trapped 
ultra-cold clouds of alkali atoms fl , various experiments 
have proved the presence of phase coherence in trapped 
condensates. The MIT group has found the inter- 
ference of two independently prepared condensates, once 
they expand and overlap after switching off the traps. 
The MIT [§, NIST |j and Munich f| experiments pro- 
vide evidence for the phase coherence of trapped conden- 
sates through the measurement of the phase coherence 
length and/or single particle correlations. 

These results support the usual picture of BEC in 3D 
gases. In equilibrium, the fluctuations of density and 
phase are important only in a narrow temperature range 
near the BEC transition temperature T c . Outside this re- 
gion, the fluctuations are suppressed and the condensate 
is phase coherent. This picture precludes the interesting 
physics of phase- fluctuating condensates, which is present 
in 2D and ID systems (see HQ and refs. therein). 

In this Letter we show that the phase coherence prop- 
erties of 3D Thomas-Fermi (TF) condensates depend on 
their shape. In very elongated 3D condensates, the axial 
phase fluctuations are found to manifest themselves even 
at temperatures far below T c . Then, as the density fluc- 
tuations are suppressed, the equilibrium state will be a 
condensate with fluctuating phase (quasicondensate) sim- 
ilar to that in ID trapped gases |7J. Decreasing T below 
a sufficiently low temperature, the 3D quasicondensate 
gradually turns into a true condensate. 

The presence and the temperature dependence of ax- 
ial phase fluctuations in sufficiently elongated 3D con- 
densates suggests a principle of thermometry for Bose- 
condensed gases with indiscernible thermal clouds. The 
idea is to extract the temperature from a measurement 



of the axial phase fluctuations, for example by measur- 
ing the single-particle correlation function. This prin- 
ciple works for quasicondensates or for any condensate 
that can be elongated adiabatically until the phase fluc- 
tuations become observable. 

So far, axial phase fluctuations have not been measured 
in experiments with cigar-shaped condensates. We dis- 
cuss the current experimental situation and suggest how 
one should select the parameters of the cloud in order to 
observe the phase-fluctuating 3D condensates. 

We first consider a 3D Bose gas in an elongated cylin- 
drical harmonic trap and analyze the behavior of the 
single-particle correlation function. The natural assump- 
tion of the existence of a true condensate at T = 
automatically comes out of these calculations. In the 
TF regime, where the mean-field (repulsive) interpar- 
ticle interaction greatly exceeds the radial (u p ) and 
axial (u) z ) trap frequencies, the density profile of the 
zero-temperature condensate has the well-known shape 
rio(p,z) = n 0m (l - p 2 /R 2 - z 2 /L 2 ), where n 0m = fi/g is 
the maximum condensate density, with /i being the chem- 
ical potential, g = Airh 2 a/m, m the atom mass, and a > 
the scattering length. Under the condition uj p 3> w z , the 
radial size of the condensate, R = (2^/ma; 2 ) 1 / 2 , is much 
smaller than the axial size L = (2fj,/muj 2 ) 1 ^ 2 . 

Fluctuations of the density and phase of the conden- 
sate, in particular at finite T, are related to elementary 
excitations of the cloud. The density fluctuations are 
dominated by the excitations with energies of the order 
of fi. The wavelength of these excitations is much smaller 
than the radial size of the condensate. Hence, the den- 
sity fluctuations have the ordinary 3D character and are 
small. Therefore, one can write the total field opera- 
tor of atoms as ip(r) = y/n (r) exp(i</>(r)), where ^(r) 
is the operator of the phase. The single-particle correla- 
tion function is then expressed through the mean square 
fluctuations of the phase (see, e.g. 0]): 
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(^(r-Mr')) = ^o(r)no(r')exp{-([^(r, r')] 2 )/2}, (1) 
with <50(r, r') = $(r) — </>(r'). The operator $(r) is given 
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by (see, e.g., §) 



(r) = [4n (r)]- 1 /2V/+(r)a, 



(2) 



where a v is the annihilation operator of the excitation 
with quantum number (s) v and energy e„, fj = u v + v v , 
and the u, v functions of the excitations are determined 
by the Bogolyubov-de Gennes equations. 

The excitations of elongated condensates can be di- 
vided into two groups: "low energy" axial excitations 
with energies e„ < hu> p , and "high energy" excitations 
with e„ > Tiujp. The latter have 3D character as their 
wavelengths are smaller than the radial size R. There- 
fore, as in ordinary 3D condensates, these excitations can 
only provide small phase fluctuations. The "low-energy" 
axial excitations have wavelengths larger than R and ex- 
hibit a pronounced ID behavior. Hence, one expects that 
these excitations give the most important contribution to 
the long-wave axial fluctuations of the phase. 

The solution of the Bogolyubov-de Gennes equations 
for the low-energy axial modes gives the spectrum tj = 
huj z y/j(j + 3)/4 |l(J, where j is a positive integer. The 
wavefunctions ff of these modes have the form 



ft (J) = 
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Note that at any z and z' the ratio of the phase correlator 
([|) to 8 2 L is a universal function oiz/L and z'/L: 



5^z,z')] 2 ) T = 8 2 L {T)f(z/L,z'/L). 



(7) 



In Fig.l we present the function f(z/L) = f(z/L, —z/L) 
calculated numerically from Eq.M). 




\z\IL 

FIG. 1. The function f(z/L). The solid curve shows the 
result of our numerical calculation, and the dotted line is 
f(z) = 2\z\/L following from Eq.(|). 



where Pj are Jacobi polynomials. Note that the con- 
tribution of the low-energy axial excitations to the phase 
operator (|^) is independent of the radial coordinate p. 

Relying on Eqs. (||) and (^), we now calculate the 
mean square axial fluctuations of the phase at distances 
\z — z'\ -C R. As in ID trapped gases 0, the vacuum 
fluctuations are small for any realistic axial size L. The 
thermal fluctuations are determined by the equation 



[S^z,z')f 
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with Nq = (8w / l5)no m R 2 L being the number of Bosc- 
condensed particles, and Nj the equilibrium occupation 
numbers for the excitations. The main contribution to 
the sum over j in Eq.(Q) comes from several lowest ex- 
citation modes, and at temperatures T 3> Tilo z we may 
put Nj = T/tj. Then, in the central part of the cloud 
(\z\, \z'\ <C L) a straightforward calculation yields 



[5Uz,z')] 2 )t = 61\z-z>\/L, 



(5) 



where the quantity 8 2 L represents the phase fluctuations 
on a distance scale \z — z'\ ~ L and is given by 

£!(T) = 32 M T/157V (^) 2 . (6) 



The phase fluctuations decrease with temperature. As 
the TF chemical potential is fi= (157Voff/7r) 2/5 (ma; 2 /8) 3 / 5 
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), Eq.(^|) can be rewritten in the form 



Si 



{T/T C ){N/N f/H 2 c 



(8) 



where T c N^-^hui is the BEC transition temperature, 
and N the total number of particles. The presence of 
the 3D BEC transition in elongated traps requires the 
inequality T c 3> Tiojp and, hence, limits the aspect ratio 
to uj p /u; z <C N. The parameter S 2 is given by 

r2 _ 32/1(^0 = N) /^\ 4/3 a 2 /5 TO i/5^2/i 5 



Except for a narrow interval of temperatures just be- 
low T c , the fraction of non-condensed atoms is small and 
Eq.(|) reduces to S 2 L = (T/T C )S 2 . Thus, the phase fluctu- 
ations can be important at large values of the parameter 
S 2 , whereas for i c 2 < 1 they are small on any distance 
scale and one has a true Bose-Einstein condensate. In 
Fig. 2 we present the quantity 5 2 for the parameters of 
various experiments with elongated condensates. 

In the Konstanz [|l3| and Hannover [Q experiments 
the ratio T/T c was smaller than 0.5. In the recent ex- 
periment JL7J the value S 2 ~ 3 has been reached, but the 
temperature was very low. Hence, the axial phase fluc- 
tuations were rather small in these experiments and they 
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were dealing with true condensates. The last statement 
also holds for the ENS experiment JjlJ where T was close 
to T c and the Bose-condensed fraction was Nq/N w 0.75. 

The single-particle correlation function is determined 
by Eq.(|l|) only if the condensate density no is much larger 
than the density of non-condensed atoms, n' . Otherwise, 
this equation should be completed by terms describing 
correlations in the thermal cloud. However, irrespective 
of the relation between no and n', Eq.(|l|) and Eqs. (|j)- 
(Q) correctly describe phase correlations in the conden- 
sate as long as the fluctuations of the condensate density 
are suppressed. This is still the case for T close to T c 
and iVo <C N, if we do not enter the region of critical 
fluctuations. Then, Eq.(|) gives 8\ = (A/A ) 3 / 5 . At the 
highest temperatures of the Bose-condensed cloud in the 
MIT sodium experiment , the condensed fraction was 
Nq/N ~ 0.1 and the phase fluctuations were still small. 
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FIG. 2. The parameter $1 for experiments with elongated 
condensates. The up and down triangles stand for <5j? in the 
sodium [11] and hydrogen [12] MIT experiments, respectively. 
Square, cross, diamond, circle, and star show 51 for the ru- 
bidium experiments at Konstanz [13], Munich [5], Hannover 
[14], ENS [15], and AMOLF [16]. 

On the contrary, for Nq/N « 0.06 in the hydrogen 
experiment p2[ , with 81 from Fig. 2 we estimate 8 2 L ~1. 
The same or even larger value of 8\ was reached in the 
Munich Rb experiment |^| where the gas temperature 
was varying in a wide interval around T c . In the Rb 
experiment at AMOLF [|l6| , the smallest observed Bose- 
condensed fraction was Nq/N 0.03, which corresponds 
to 8\ « 5. However, axial fluctuations of the phase have 
not been measured in these experiments. 

We will focus our attention on the case where N 
N and the presence of the axial phase fluctuations is 
governed by the parameter 5%. For 81 3> 1, the nature 
of the Bose-condensed state depends on temperature. In 
this case we can introduce a characteristic temperature 



at which the quantity 8\ « 1 (for Ao ~ N). In the 
temperature interval < T < T c , the phase fluctuates 
on a distance scale smaller than L. Thus, as the density 
fluctuations are suppressed, the Bose-condensed state is 
a condensate with fluctuating phase or quasicondensate. 
The expression for the radius of phase fluctuations (phase 
coherence length) follows from Eq.(||) and is given by 



14, « L(T /T). 



(11) 



T = 15(?k^) 2 A/32u 



(10) 



The phase coherence length 1$ greatly exceeds the cor- 
relation length l c = H/^/mfi. Eqs. (|ll]) and (10) give 
the ratio l^/lc w {T c /T)(T c /Hlu p ) 2 > 1. Therefore, the 
quasicondensate has the same density profile and local 
correlation properties as the true condensate. However, 
the phase coherence properties of quasicondensates will 
be drastically different (see below). 

The decrease of temperature to well below makes 
the phase fluctuations small (<5| <C 1) and continuously 
transforms the quasicondensate into a true condensate. 

It is interesting to compare the described behavior of 
the interacting gas for Si ^> 1, with the two-step BEC 
predicted for the ideal Bose gas in elongated traps [fT8| . 
In both cases, at T c the particles Bose-condense in the 
ground state of the radial motion. However, the ideal gas 
remains non-condensed (thermal) in the axial direction 
for T > T 1D = Nhu> z /\n2N (assuming T 1D < T c ), and 
there is a sharp cross-over to the axial BEC regime at 
T Tid- The interacting Bose gas below T c forms the 
3D TF (non-fluctuating) density profile, and the spatial 
correlations become non-classical in all directions. For 
81 1, the axial phase fluctuations at T ~ T c are still 
large, and one has a quasicondensate which continuously 
transforms into a true condensate at T below T^. Note 
that is quite different from T\u of the ideal gas. 

Let us now demonstrate that 3D elongated quasi- 
condensates can be achieved for realistic parameters of 
trapped gases. As found above, the existence of a qua- 
sicondensate requires large values of the parameter 81 
given by Eq. @ . Most important is the dependence of 81 
on the aspect ratio of the cloud w p /w z , whereas the de- 
pendence on the number of atoms and on the scattering 
length is comparatively weak. Fig. 3 shows T c /T<^ = 81, 
/i/T^, and the temperature as functions of uj p /uj z for 
rubidium condensates at N = 10 5 and lo p = 500 Hz. 
Comparing the results for Si in Fig. 3 with the data in 
Fig. 2, we see that 3D quasicondensates can be obtained 
by transforming the presently achieved BEC's to more 
elongated geometries corresponding to uj p /uj z > 50. 

One can distinguish between quasicondensates and 
true BEC's in various types of experiments. By using the 
Bragg spectroscopy method developed at MIT one can 
measure the momentum distribution of particles in the 
trapped gas and extract the coherence length l^, || . The 
use of two (axially) counter-propagating laser beams to 
absorb a photon from one beam and emit it into the other 
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one, results in axial momentum transfer to the atoms 
which have momenta at Doppler shifted resonance with 
the beams. These atoms form a small cloud which will 
axially separate from the rest of the sample provided the 
mean free path greatly exceeds the axial size L. The 
latter condition can be assured by applying the Bragg 
excitation after abruptly switching off the radial confine- 
ment of the trap. The axial momentum distribution is 
then conserved if the dynamic evolution of the cloud does 
not induce axial velocities. According to the scaling ap- 
proach fll9| ], this is the case for the axial frequency de- 
creasing as Wg(t) — lj z (0)[1 + LUpt 2 ]^ 1 / 2 . 
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FIG. 3. The ratios T c /T^ = 8 2 c and fi/T^ in (a) and the 
temperature in (b), versus the aspect ratio lu p /lo z for 
trapped Rb condensates with N = 10 5 and lj p = 500 Hz. 

In "juggling" experiments described in J7J and similar 
to those at NIST and Munich , one can directly mea- 
sure the single-particle correlation function. The latter 
is obtained by repeatedly ejecting small clouds of atoms 
from the parts z and z' of the sample and averaging the 
pattern of interference between them in the detection re- 
gion over a large set of measurements. As follows from 
Eqs. (||) and (Q), for z' — —z the correlation function de- 
pends on temperature as exp {— 5 2 L {T) f (z / L) / 2} , where 
f(z/L) is given in Fig.l. 

The phase fluctuations are very sensitive to tempera- 
ture. From Fig. 3 we see that one can have T^/T c < 0.1, 
and the phase fluctuations are still significant at T < 
where only a tiny indiscernible thermal cloud is present. 



This suggests a principle for thermometry of 3D Bose- 
condensed gases with indiscernible thermal clouds. If the 
sample is not an elongated quasicondensate by itself, it is 
first transformed to this state by adiabatically increasing 
the aspect ratio lu p /uj z . This does not change the ratio 
T/T c as long as the condensate remains in the 3D TF 
regime. Second, the phase coherence length 1$ or the 
single-particle correlation function are measured. These 
quantities depend on temperature if the latter is of the 
order of or larger. One thus can measure the ratio 
T/T c for the initial cloud, which is as small as the ratio 
T^/Tc for the elongated cloud. 

We believe that the studies of phase coherence in elon- 
gated condensates will reveal many new interesting phe- 
nomena. The measurement of phase correlators will al- 
low one to study the evolution of phase coherence in the 
course of the formation of a condensate out of a non- 
equilibrium thermal cloud. 

When completing the paper we were informed that 3D 
quasicondensates were observed in Hannover [f20f . 
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